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Abstract — We present novel exact expressions and accurate 
closed-form approximations for the level crossing rate (LCR) 
and the average fade duration (AFD) of the double Nakagami- 
m random process. These results are then used to study the 
second order statistics of multiple input multiple output (MIMO) 
keyhole fading channels with space-time block coding. Numerical 
and computer simulation examples validate the accuracy of the 
presented mathematical analysis and show the tightness of the 
proposed approximations. 

Index Terms — Level crossing rate (LCR), Average fade dura- 
tion (AFD), keyhole MIMO fading channels, Nakagami-m fading, 
multiplicative fading 



I. Introduction 

Recently, special attention has been given to the so-called 
"multiplicative" fading models. The double Rayleigh (i.e., 
Rayleigh*Rayleigh) channel fading model has been found to 
be suitable when both transmitter and receiver are moving 
(Qj. Moreover, it has also been recently used for keyhole 
channel modeling of multiple-input multiple-output (MIMO) 
systems ||2-||3l. Its extension, the double Nakagami-m (i.e., 
Nakagami-m*Nakagami-m) fading model, has been consid- 
ered in |4|, where the fading between each pair of transmit 
and receive antennas in presence of the "keyhole" is charac- 
terized as Nakagami-TO fading. However, all the above works 
describe and utilize only the first order statistical properties of 
these "multiplicative" fading models, such as the outage and 
the error probabilities. But, knowledge of the second order 
statistics for above fading models are equally important, and 
are applicable, for example, in modeling and design of the 
multihop communications systems |5|. 

In this letter, we focus on the second order statistics 
of the double Nakagami-m random process, for which we 
determine exact and approximate analytical solutions for its 
level crossing rate (LCR) and average fade duration (AFD). 
Then we apply these results to study the second order statistics 
of the keyhole channels applicable to MIMO systems with 
space-time block coding (STBC), operating in specific rich- 
scattering environments. Note that although this work assumes 
independence among the channels, similar analysis can be used 
to derive LCR and AFD in correlated keyhole channels |9|. 
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II. On the second order statistics of the double 

NAKAGAMI-m RANDOM PROCESS 

Let the double Nakagami-m random process be defined as 



Z{t) = X{t)Y{t) , 



(1) 



where X{t) and Y{t) are a pair of independent Nakagami-m 
distributed RVs with probability distribution functions (PDFs) 
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where fix — E[X^], Qy — E[Y^], and mx and my are the 
fading severity parameters, where E[-] means expectation. 

If X{t) and Y{t) are signal envelopes in some scattering 
radio channel exposed to the Doppler effect due to stations' 
relative mobility, then X{t) and Y{t) are time-correlated 
random processes. Considering a fixed-to-mobile channel, 
each scattered component of X{t) and Y{t) has some result- 
ing Doppler spectra with maximum Doppler frequency shift 
fmx and f„iy, respectively. It was shown in |6| that, under 
such conditions, the envelopes time derivatives X and Y are 
independent from their respective envelopes, while following 
zero-mean Gaussian PDFs with respective variances 



= (TTfrnx) i^x/mx, 
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A. Second order statistics 

The LCR of Z at threshold z is defined as the rate at which 
the random process crosses level z in the negative direction. 
To extract LCR, we need to determine the joint PDF of Z and 
Z, f^^{z,z), and apply the Rice's formula 



Nziz)^ / zfz^{z,z)dz. (5) 
Jo 

The above expression can be rewritten as 

Nz{z)=J (j zf^^^xi^\z,x)dzjfz\x{Ax)fx{x)dx 

° ° . (6) 

where fz\zxi'^ ' ') conditional PDF of Z conditioned 

on Z and X. This conditional PDF can be determined by 
finding the time derivative of both sides of ([U, 
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from which it is easily seen that, for fixed Z ~ z and 
X = X, the time derivative Z is a zero-mean Gaussian RV 

.2 2_2 /_2 I ^2^2 



integral in ^ can be solved as 



The critical point of f{x) is determining as the value of x for 
which df /dx = 0, i.e., 



with variance ct;; 
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The conditional PDF of Z for some fixed X — x, fz\xiz\x), 
is determined by simple transformation of RVs, fz\x{z\x) = 
fY{z/x)/x. Substituting dHJ into (|6]l, after some algebraic 
manipulations, we obtain the exact solution for the LCR 
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Using (fTsll-dTsll. the approximate closed-form solutions for the 
LCR and the AFD are respectively obtained as 
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The above integral can be evaluated numerically with de- 
sired accuracy (e.g. by using some common software such 
as Mathematica). Alternatively, one can apply the Laplace 
approximation to obtain a highly accurate closed-form solution 
of (|9]l - as presented in the following subsection. 

The AFD of Z at threshold z is defined as the average time 
that the double Nakagami-m random process remains below 
level z after crossing that level in the downward direction. 



Tz{z) 



Fz{z) 



(10) 



Nz{z) ' 

where Fz{z) denotes the CDF of Z, which was derived only 
recently in closed-form for A^*Nakagami random process Q. 
For the double Nakagami random process, it attains the form 
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where r(-) and G[-] are gamma and Meijer's G functions. 

B. Laplace approximation 

Using IHl, the Laplace type integral can be approximated as 



g{x) e-^f'^^'Ux : 
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(12) 



when the real valued parameter A is very large (i.e., A 
oo). In (fTzl l. f{x) and g{x) are real-valued functions of x 
and xo is the point at which f{x) has an absolute minimum 
(known as the interior critical point of f{x)). Note, that f"{x) 
denotes the second derivative of f{x) with respect to x. It was 
observed that above approximation is very accurate even for 
small values of A |8|. Comparing (fT2l i and (|9]l, these functions 
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whereas the second derivative of the former is / (x) — 

2TOx/f^x + 6(myz^)/(riya;'*) + 2(mx — TOy)/a:^ and A = 1. 
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Although substitution of f{xo), f (x^) and g(xo) into (fTSI l 
and ( [TtI i is omitted for brevity, we emphasize that the threshold 
z appears only as the ratio z'^ / {{ilx /mx){^Y /mY))- 

III. MIMO STBC Communication over Keyhole 
Fading Channels 

Potentials of MIMO communications systems are not al- 
ways achievable even for a fully uncorrelated transmit and 
receive channels, which is attributed to the rank deficiency of 
the MIMO channels known as the keyhole or pinhole effect 
f2]. The existence of the keyhole MIMO channels has been 
proposed and demonstrated through physical examples, where, 
although spatially uncorrelated, these channels still have a 
single degree of freedom |2|-|3|. Under the keyhole effect, the 
entries of the channel matrix, H, follow statistics described as 
a product of two independent single-path gains. 

A. The MIMO keyhole channel model 

From [4], the complex path gain of baseband equivalent 
signal transmitted over the channel between the i-th transmit 
and the j-th receive antenna at arbitrary moment t is expressed 

as 1 < i < 1 < j < 



hij{t) 



(18) 



where {ai(i)e-'"^''^*^}^j^ are the complex path gains intro- 
duced by the rich-scattered channel from the i-th transmitting 



antenna to the "keyhole", and {/3j(t)e^'''^^*''}^^ are the com- 
plex path gains introduced by the rich-scattered channel from 
the "keyhole" to the j-th receiving antenna. Phases {4'i{t)}f^-^ 
and {ipj{t)}^^^ are independent and uniformly distributed 



over [0, 27r). The amplitudes {ai{t)}fLi and {I3j{t)}"^^ are 
i.i.d. Nakagami-TO RVs. The fading severity parameters of 
ai{t) are equal to mx, whereas fir — E[af] for all i. Simi- 
larly, the fading severity parameters of f3j (t) are equal to to/?, 
(^^-^ whereas — E[Pj] for all j. Assuming mobility of both 
the transmitter and the receiver with respect to the "keyhole", 
all channel gains are time-correlated random processes with 
maximum Doppler shifts /q, . — fa and fp. — fp, respectively. 



I N 



3 



1.2 


Approximalion 

o UGsimulalion 








M = 7.N = 




0.8 

0.6 


















2, N = 2?/ 










0.4 








M = 3, N = 

W 








0.2 














-°, 














5 -1 









10 15 20 



Normalized threshold (dB) 



Fig. 1 . Normalized LCR for various number of transmit and receive antennas 

Under such conditions, the time derivatives and 0j are 
independent from and /3j, respectively, and both follow 
zero-mean Gaussian PDFs with variances given by (|4|i, cr? = 

{■KfcfnT/mT, 1 < i < M and ct^ ^ {■KfpfQ.p^/niR, 
l<j <N. 

B. Orthogonal space-time block coding and decoding 

The orthogonal spacetime block encoding and decoding 
(signal combining) transform a MIMO fading channel into 
an equivalent single-input-single-output (SISO) fading channel 
with a path gain of the squared Frobenius norm of the MIMO 
channel matrix H(t) = [hij{t)]j^j^j^ 14], 

M N / M \ ( ^ \ 

\mmi = EE I'^'.wi' = E^'w E/^lw 

(19) 

at arbitrary moment t. After space-time block decoding, the 
instantaneous output signal-to-noise ratio (SNR) per symbol 
is given by 

iit) = ^\mt)\?F, (20) 

where 7 = Es/Nq is the average SNR per receive antenna, 
and R is the rate of the STBC. 

C. Second order statistics of output SNR 

We introduce the auxiUary random process Z{t) defined by 

Z{t)^^\m)\\l^X{t)Y{t). (21) 

where Xit) = ySXW and Y{t) = 
again Nakagami-m distributed with PDFs given by (|2]l and 
(|3]l, respectively, with mx = MniT, — MUt, my = 
Nmji and ily = NVlj^. The time derivatives X and Y are 
independent from X and Y, respectively, and both follow the 
zero-mean Gaussian PDF with variances given by (|4]i, cr^ = 
^ii = {-Kfaf^T/rriT and a\ = aj = [TrfpfVLBj-mR. 

Hence, the random process Z{t], defined by (I2TI 1. is a 
double Nakagami-TO process for which we can apply the 
analytical framework of Section II to determine its exact and 
approximate LCR and AFD by using (|9]l, ([TTli, ([T6]l and ( [TtT i. 
With above in mind, the LCR and the AOEU of instantaneous 

'instead of tlie term "average fade duration (AFD)", tlie term "average 
outage duration (AOD)" is used here. 
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Fig. 2. Normalized AOD for various number of transmit and receive antennas 

output SNR, given by ( |20] |. are respectively determined as 

N-,[l) - Nz{ VjMR/j ), (22) 
r^(7) = TzWjMR/j). (23) 

IV. Numerical Results 

We present several numerical examples for the LCR and the 
AFD of the STBC MIMO communications system operating 
over a keyhole fading channel. The mobile transmitter and 
the mobile receiver are assumed to introduce same maximum 
Doppler shifts due to same relative speeds with respect to the 
"keyhole", yielding = = f„,. 

Figs. 1 and 2 depict the normalized LCR {N^/fm) and 
normalized AFD (T^ /„) of the instantaneous output SNR vs. 
normalized SNR threshold. The normalized SNR threshold (x- 
axis) is calculated as lQ\og[j M R/ {^{Ht /mT){i^R/rnii))]. 
The results are obtained for three different pairs of number of 
transmit and receive antennas {M, TV), appearing as curve pa- 
rameters. For each pair (A/, N), the three comparative curves 
on both figures indicate excellent match between the exact and 
the approximate solutions for the two statistical parameters, 
both of which are validated by Monte Carlo simulations. 

References 

[1] J. B. Andersen, "Statistical distributions in mobile communications using 
multiple scattering," Pmc. Gen. Assem. Int. Union of Radio Sci. , Maas- 
tricht, The Netherlands, Aug. 2002. 

[2] D. Gesbert, H. Bolcskei, D. A. Gore, and A. J. Paulraj, "Mimo wireless 
channels: Capacity and performance prediction," Pmc. GLOBECOM 2000, 
vol. 2, pp. 1083-1088 

[3] D. Chizhik, G. J. Foschini, M. J. Gans, and R. A. Valenzuela, "Keyholes, 
correlations, and capacities of multielement transmit and receive antennas," 
IEEE Trans. Wireless Commun., vol. 1, no. 2, pp. 361-368, Apr 2002 

[4] H. Shin and J. H. Lee, "Performance analysis of space-time block codes 
over keyhole Nakagami-m fading channels," IEEE Trans. Veh. Techno!., 
vol. 53, no. 2, pp. 351-362, Mar. 2004 

[5] Z. Hadzi-Velkov, N. Zlatanov, G. K. Karagiannidis, "On the second 
order statistics of the multihop Rayleigh fading channel," accepted for 
publication in the IEEE Transactions on Communications, 2008 

[6] M. D. Yacoub, J. E.V. Bautista, and L. G. de Rezende Guedes, "On higher 
order statistics of the Nakagami-m distribution," IEEE Trans. Commun., 
vol. 48, No. 3, pp. 790-794, May 1999 

[7] G. K. Karagiannidis, N. C. Sagias, and P. T. Mathiopoulos, 
"7V*Nakagami: a novel stochastic model for cascaded fading fhannel," 
IEEE Trans. Commun., vol. 55, No. 8, pp. 1453-1458, Aug. 2007 

[8] R. Wong, Asymptotic Approximations of Integrals, SIAM: Society for 
Industrial and Applied Mathematics, New edition, 2001. 

[9] G. Levin, S. Loyka, "On the Outage Capacity Distribution of Conelated 
Keyhole MIMO Channels", IEEE Trans, on Inform. Theory, vol. 54, no. 7, 
pp. 3232-3245. July 2008. 



